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Abstract 

The complete tree level cross-section for the process e + e~ — > e~V e ud 
is computed using the GRACE system, a program package for auto- 
matic amplitude calculation. Special attention is brought to the gauge 
violation problem induced by the finite width of the VF-boson. The 
preserved gauge scheme is introduced and an event generator including 
double-resonant, single-resonant and non-resonant diagrams with no 
need for a cut on the electron polar angle is built. Below threshold, the 
single W and non-resonant diagrams give a substantial contribution 
to the total cross-section, at linear collider energies, the cross-section, 
for small electron polar angles, is simply dominated by these contribu- 
tions. Since no cut needs to be applied to the electron, the generator 
can be used to estimate background for searches involving jets and 
missing energy. A mono jet event rate estimation based on this pro- 
cess at LEP-I energy is discussed. 



To be submitted to Physics Letters B. 
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1 Introduction 



Precise event generators for processes involving W, Z, 7 and leading to four 
fermions in the final state are necessary for a good understanding of the 
three-boson coupling at LEP and at future linear collider CM. energies. 
Four-fermion final state can be produced by double (heavy boson) resonant 
diagrams, single resonant diagrams or non-resonant diagrams. Single W pro- 
cesses include e + e~ — > l~v{W + where I = e, fi. The I = e case differs from 
the I = fi one by the existence of the diagrams where 7, Z° or even W are 
exchanged in the t-channel. Most of these diagrams have a 7-7 like behavior 
such as a strong forward electron peak, the difference lying in the propaga- 
tor mass and width. These diagrams will be called the t-channel diagrams 
hereafter in this paper. Since the VT-boson is unstable and decays, final 
state like e + e~ — > e~v e ud will be actually observed. The diagrams involved 
in e + e~ — > e~v e ud can be grouped into the s-channel (Fig.l) and t-channel 
classes (Fig. 2), each group forms a gauge invariant set. In regard to the W- 
resonance, they can also be characterized as double-resonant, single-resonant, 
and non-resonant diagrams. The single and non-resonant diagrams play an 
important role on several studies, including the W mass determination, the 
anomalous couplings limits and the background estimates in the search of 
new particles. 

Recently event generators for four-fermion (+ photons) processes have 
been developed ||. However, in order to cope with the complexity of the 
calculations, some approximations are always introduced. On the contrary, 
the GRACE package can automatically produce the complete set of the tree 
level diagrams involved in the process while taking into account fermion 
masses. 

In a previous work H, a complete calculation of two typical processes 
of four-fermion final states in e + e~ collisions, e + e~ — > e~v e ud and e + e _ — > 
udud, has been presented for the first time. We have pointed out that the 
single and non-resonant diagrams play an important role below the iy-pair 
threshold. Their contribution, for instance, reaches 27% for the udud process 
for y/s = 150 GeV and gives a non-negligible effect even at higher energy 
(4.4% for e~v e ud at yfs = 190 GeV). Below the threshold only the off-shell 
double-resonant amplitudes compete with the single- and non-resonant ones, 
so that the relative contribution of the latter becomes quite large. Besides 
a set of experimental loose constraints, a cut on the polar angle of the final 
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electron with respect to the initial electron direction, 6 e , was applied to avoid 
a significant gauge violation appearing in the subset of j-W diagrams. 

In this paper, we analyze the effect of the gauge violating term due to 
the finite width of the W^-boson. The so-called preserved gauge scheme is 
introduced to overcome this problem. The total e + e~ — > e~v e ud cross-section 
from below the P^-pair threshold up to linear collider energy, with no cut on 
the final electron, is presented. As no constraint need to be applied on the 
electron kinematic, this generator can be used to estimate the background 
to signal requesting jets and missing energy like the search for new particles. 
As an example, the contribution of this process for mono-jet production at 
LEP-I is estimated in the last section of this paper. 

2 The GRACE system 

The GRACE system |§ has been developed to perform the very lengthy 
computations involved in the study of high energy reactions. The GRACE 
package is a complete set of tools for computing tree level processes. All the 
usual steps occurring in a given computation are covered: from the process 
specification to the event generator. It is composed of three components: the 
diagram generator, the matrix element builder using helicity amplitudes from 
the CHANEL Q library and the multi-dimensional phase space integration 
package BASES || associated with the event generator SPRING 0. 

Fermion masses are properly taken into account in the helicity ampli- 
tudes. The boson width is introduced into the gauge propagator when the 
denominator may vanish for positive squared momentum transfer. A gauge 
invariance checking program is automatically built by the system. 

The results presented hereafter have been obtained using the following 
set of parameters: 



M z 


= 91.1 GeV 


T z 


= 2.534 GeV 


a 


= 1/137 


sin 2 9 W 


= 1 - (Af^/Aff) 


M w 


= 80 GeV 


m u 


= m d = 0.1 GeV. 
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The W width is taken from the Particle Data Group Table: IV = 2.25 
GeV. The gauge boson (W, Z) widths are assumed to be constant in the 
calculation. Furthermore some realistic experimental cuts have been intro- 
duced: 



0° < B e - < 180° case - a 
8° < B e - < 172° case - b 



8° < 9 u j < 172° , E u j > 1 GeV 



where 6 e - is the final state electron polar angle, measured from the incident 
e~ beam, B u j are the similar angle for the u and d quark and E u g are the 
energies of final u and d quarks. 

The gauge invariance of the amplitude without particle width is checked 
numerically by a random selection of the boson gauge parameters at several 
points of the phase space. The errors are within the precision of the numerical 
calculation (typically less than C(10~ 12 ) in double precision). 



3 Effect of the gauge violating term 

The violation of the gauge invariance at the tree level is due to the intro- 
duction of the IV-boson finite width. The first four diagrams in Fig. 2, the 
so-called j-W diagrams, give the dominant contribution to the t-channel 
amplitude, the e~ being scattered in the forward direction. Since a large 
cancellation occurs among 7- W diagrams , the gauge violating terms lead 
to a strongly divergent cross-section at small electron polar angle. It blows 
up by about six orders of magnitude at B e ~ when the width is introduced 
directly in the propagator (Fig.3, dashed line). 

To see how the effect of the gauge violating term arises, we examine the 
total amplitude of the j-W diagrams. Since there is only one W propagator 
whose four-momentum transfer squared is positive, the amplitude without 
the H^-width can be written as: 

M = -pV„ (1) 

^ = u(p'hMp)[9vu + (f - tykpK/k 2 ], (3) 
= u(p'hMp)> ( 4 ) 
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where p^(p^) is the four-momentum of the initial (final) electron, = Pa—p' a i 
the momentum of the virtual photon, g M , the J^-propagator four-momentum 
transfer (q 2 > 0) and r^n^), the single-resonant (non-resonant) diagram 
amplitudes. If the width, T w , is introduced, one obtains: 

T ^ T '» = q 2 -M 2 w l iM w T w + (5) 
d a + iM w Y w n a ^ 



q 2 -M 2 v + iM w Y w ' 
dn = r^ + (q 2 - M^)n n . (7) 

The square of the electron current, after averaged over spin states, is: 

= ^I/mK = 2[pX + PuP'„ + -j9^\- ( 8 ) 

spin 

Then the squared amplitude is given by: 

IA/fl 2 2e 2 (L^d^dl) + M^Y 2 v iL ViV n v ji* v ) . . 

' ' " (k 2 ) 2 {q 2 - M^) 2 + M 2 W T 2 W ' W 

(note that and n M are real numbers at the tree level). The first term in 
numerator is the gauge invariant part; when k 2 — > 0, it behave as L nu d^dl — > 
0(k 2 ). The second term is the gauge violating one as L^n^n* — > 0(1) when 
k 2 — > 0. By integrating over /c 2 , the former gives the well-known log(s/m 2 ) 
dependence of the total cross-section. The latter, the gauge violating term, 
however, dominates the total cross-section after integration, as it does not 
compensate the photon propagator oc l/(k 2 ) 2 . It is clear that this term is 
responsible for the divergence of the cross-section in the small angle region 
(Fig. 3, dashed line). It should be emphasized that this behavior does not 
depend on the gauge parameter of the W propagator. 

The introduction of the W width, yet preserving the gauge invariance of 
the amplitude, can be achieved by applying the following method. 

The current l n is replaced by the momentum k^, gauge invariance implies: 

k ^ = if 2 + k - n ' ( 10 ) 

q — ivi w 
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Then is a gauge invariant quantity. T M can be expressed in term of as: 

If the particle width is introduced at this stage, the amplitude can be cast 
into the form: 

T — > T" = ^ fis) 

" " q 2 -M^ + iM w r w {i6> 

This amplitude is apparently gauge invariant as the divergent term discussed 
previously has disappeared || . The total cross-sections based on this ampli- 
tude do not diverge even in the small angle region as shown by the solid line 
in Fig. 3. This scheme, the preserved gauge scheme, can be interpreted as fol- 
lows; the tree-level amplitude is gauge invariant but divergent at q 2 = M^- 
To avoid this divergence, an imaginary part of a higher-order is introduced 
into the propagator denominator as a particle width. However this order mix- 
ing causes the violation of the gauge invariance at the tree-level and gives 
rise to the cross-section divergence. To preserve the gauge invariance, an ad- 
ditional term is added to the amplitude on the analogy of the counterterms, 
it gives the higher order correction to the leading term. 

In order to check the validity of the method, one can approach this result 
by applying a method which minimize the gauge violating term. Let us 
remind a usual trick used for the tensor L^ v . Let's assume, first, that T M is 
gauge invariant. One can replace L^ u by: 

L p.v -> L'^ = Ap^p u + k 2 g^ u . (14) 

In this equation, the first term is responsible for the blow-up of the cross 
section. Because of the gauge invariance, one can further replace the vector 
Pii by: 

Vn -> Pn = P M - (Po/k )k^ ( 15 ) 

where po and k are the 0-th components of the four- momenta p^ and k^, 
respectively. By substituting p^ in Eq.(15) and dropping k^, one gets 

J> - L % = ±P»P» + fcV- (16) 

It is known that a product P ■ A of P with an arbitrary vector A can be 
expressed by a sum of terms proportional to either m 2 ., 1 — cos# e or sin# e . 
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Hence in the region 9 e ps 0, both P ■ d and P ■ n behave like k 2 , because 
1 — cos 9 e vanishes almost like k 2 . 

Then, one possible way to get rid of the large gauge violation is to use the 
current L" , instead of the original one in Eq.(9). The new amplitude 
squared becomes: 

|yu|2 2e 2 (L'^dl) + M 2 w T 2 w (L%n,nl) 

' ' ~" (k 2 ) 2 {q 2 - M 2 W ) 2 + M 2 W T 2 W ■ 1 l} 

The gauge violating term L'^ v n^n u behaves as 0(k 2 ) for vanishing k 2 , and it 
should not induce violent divergence. Since the non-resonant diagrams have 
a lU-boson space-like propagator, one may assume roughly that (9(r M ) ~ 
C(n M M^). Hence the gauge violating term is proportional to (T w /M w ) 2 ~ 
10~ 3 . This is small compared to the gauge invariant term (L^d^d*) and can 
be neglected. 

Using this last method, the cross-section turns out to be 8.63xlCT 2 pb 
(case-a at y/s = 180 GeV) as seen in Fig. 3 (indicated by the arrow). This is 
only 5% larger than the cross-section obtained by the preserved gauge scheme 
(8.24xl0- 2 pb). 



4 Threshold behavior 

In the following, the computations are performed using the preserved gauge 
scheme. The contribution from the t-channel diagrams strongly depends on 
6 e cut as shown in Fig. 4. If we require the electron polar angle to be in the 
range 172° > 9 e > 8°(140 mrad) case-b, the double-resonant diagrams are 
dominant (96% at y/s = 180 GeV) but their contributions are reduced down 
to 83% if the electron angular cut is set to a vanishing value (case-a). 

The threshold behavior of e + e~ — > e~v e ud total cross-section is presented 
in Fig. 5 together with one high energy point (y/s = 500 GeV). With no 
electron angular cut, case-a, one observes a large increase in the cross-section 
due to the t-channel diagrams, although when such a cut is applied, case-b, 
the effect remains substantial only above threshold. The relative contribution 
of the double-resonant diagrams is shown in Fig.6. At y/s = 500 GeV, the 
share of the double-resonant diagrams is only 17% without electron tagging 
(case-a), while it amounts to 55% when the electron is tagged (case-b). 
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For the W mass measurement, it should be noted that using the threshold 
scanning method, a precise electron tagging is necessary to deal with the t- 
channel contribution. When the direct reconstruction of 2-jets is used, the 
invariant mass distribution of a quark pair gets a low mass tail from the 
non-resonant diagrams. The effect of the electron angular cut is to suppress 
the contribution from those diagrams, as shown in Fig. 7. 

In the studies of the anomalous couplings of the gauge bosons, the pro- 
cess e + e~ — > eVeW is considered. However the process to be measured is 
e + e~ — > ev e qq [f|J which includes the effect of the non-resonant diagrams. 
The contributions of i-channel diagrams (Fig. 2) have been computed for both 
processes e + e~ — > eu e W and e + e~ — > eu e qq, with no cut, a 14% discrepancy 
is found in the total cross-section at yfs = 180 GeV, which comes from the 
small angle region of the final electron. 

alee — > euW) * Br(W — > ud) 
1 — = 0.14 

o[ee — > evud) 

This effect should definitely be taken into account in these studies. 

The preserved gauge scheme can also be applied to different processes, 
like those appearing at electron-electron collider. For example, the process 
e~e~ — > e~W~v e |J involves only the i-channel diagrams at the tree level. 

A comparison with the four-fermion Monte-Carlo EXCALIBUR |]TD has 



been performed for the e + e~ — ► e~v e ud channel as an independent test of 
our work. However such a comparison make sense only when a cut is applied 
on the electron polar angle as EXCALIBUR does not treat the forward 
divergence. Under this condition, perfect agreement is obtained with our 
generator without the preserved gauge scheme. Identical divergent behavior 
is found when the cut is decreased down to almost zero. The net result of our 
study is the extension of the domain of validity of the generator for massive 
fermions and for vanishing electron polar angle. 



5 Monojet event rate at LEP 

The ALEPH collaboration at LEP has reported the observation of three 
monojet events, one leptonic and two hadronic, from a data sample of 82 
P^ 1 EH recorded at and close to the Z° peak. Although the event rate is 
consistent with the expectation from e + e~ —>■ ^*vv with 7* — > //, the large 
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observed mono jet mass and transverse momentum are quite unlikely in this 
process. Even when diagrams from Z or W exchanges are taken into account, 
the probability of their occurrence is about 5%. 

For the W exchange process, e + e~ — > euff', only four types of diagrams 
are taken into account, double-resonant and single-resonant diagrams in s- 
channel and single- resonant diagrams in t-channel. Non-resonant diagrams 
were not included as no generator was available at that time. These diagrams 
are expected to give a large contribution to the total cross-section at energies 
below the W^-pair threshold, moreover the mass and transverse momentum 
distribution are expected to peak at higher values. 

The expected number of events from e + e~ — *> e ±<y u\qq' and the probabil- 
ity to produce the observed jet masses and transverse momenta have been 
estimated by applying the ALEPH cuts at the generator level. 

1. 9 e < 25 mrad. 

2. Polar angle of qq' system < 25.8°. 

3. Neither q nor q' goes backward of qq' system. 

4. > 1.3 GeV. 

The total cross-section of the e + e~ — > e^v^qq' processes (qq' = ud and cs) 
at y/s = 91.1 GeV with the above cuts is 7.2 x 10~ 3 pb, which corresponds 
to 0.59 events for a integrated luminosity of 82 pb~ l . The expected event 
distribution of qq' invariant masses and transverse momenta are shown in 
Figs. 8 and 9, respectively with the same binning as Fig. 3 in ref. ||ll|| . The 



observed jet masses and transverse momenta of the hadronic monojet (M = 
3.2, 5.3 GeV and pr = 6.6, 18.5 GeV) are distributed around the peak 



of the distributions. The probability defined in ref. |TTJ is calculated to be 
0.95. The observed events are compatible with the expectation obtained from 
e + e~ — > e ±<y u\qq' processes. However detailed calculations with hadronization 
and detector simulation are needed to give a more precise estimation. 



6 Summary and Conclusions 

The complete tree level cross-section of the process e + e~ — > e~v e ud has been 
computed using the GRACE package. It is shown that the naive method of 
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introducing the W width induces an unphysical divergent total cross-section 
due to the gauge violation among the t-channel tree-level diagrams. This 
effect is sizable when the electron scattering angle is less than a few degrees. 
Two ad hoc methods are discussed, one curing definitely the violent blow 
up of the amplitude (the preserved gauge scheme) and the other strongly 
reducing this violation. 

The cross-sections as a function of the CM energies without any cut on the 
final electron are presented. We conclude that the approximation based only 
on a part of the diagrams (the double-resonant diagrams and/or the single- 
resonant '-f-W diagrams) is not sufficient to give precise predictions neither 
at threshold energy nor at high energy. For the W mass determination 
based on the W pair threshold scanning, the complete calculation should 
be used and a precise measurement of the electron should be made. The 
single-resonant diagrams can be used for the W mass measurement by the 
2 jet reconstruction method, although the kinematic constraints are much 
weaker. Observed jet masses and transverse momenta of hadronic monojet 
seen by the ALEPH collaboration at LEP are not in contradiction with the 
expectation from e + e~ — > e^v\qq' processes when naive experimental cuts 
are applied at parton level. 
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Figure 1: The s-channel diagrams of the e + e — > e v e ud process in unitary- 
gauge. The first three diagrams in the first row are double-resonant diagrams. 
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Figure 2: The t-channel diagrams of e + e~ — > e~~v e ud process in unitary gauge. 
The first and second columns show the single-resonant diagrams and the rest 
shows the non-resonant diagrams. Diagrams in first row(7-l¥ processes) give 
the dominant contribution among t-channel diagrams. 
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Figure 3: The cross-section vs. electron cut angle for the sub-set of ^y-W 
diagrams only at = 180 GeV. The left half of the figure is a magnified 
view of the small angle region. The dashed line shows a result with a naive 
Breit-Wigner form for the W^-propagator and the solid line corresponds to 
the introduction of the width in a gauge-invariant way using the so-called 
preserved gauge scheme as explained in the text. The arrow shows the cross- 
section using the second method described in the section 3. 
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Electron cut angle (mrad) 

Figure 4: The cross-section vs. the e~ cut angle at \/s — 180 GeV. The solid 
line represents the contribution from all diagrams, the dashed line from the 
double-resonant diagrams, and the dotted line from the ■y-W diagrams. 
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VS (GeV) 

Figure 5: The threshold behavior of the total cross-section of the e + e" — > 
e~v e ud process. Cuts described in section 2 are applied on u and d quarks, a) 
no cut on the final electron (case-a), and b) 9 e > 8° (case-b). The solid line 
shows the cross-section from all diagrams, the dashed line from the double- 
resonant diagrams, and the dotted line from the t-channel diagrams. Results 
at y/s — 500 GeV are shown by solid arrows (all diagrams) and dashed arrows 
(double- resonant diagrams). 
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Figure 6: The relative contribution of the double-resonant diagrams to the 
cross-sections at y/s — 180 GeV. In a), no cut on the out-going electron but 
only on u and d (case-a) and in b) an additional cut, 9 e > 8° (case-b), is 
imposed. Results at y/s — 500 GeV are shown by arrows. 
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Figure 7: The invariant mass distribution of u and d calculated at ^/ s = 180 
GeV. The solid line (case-a), and the dotted one with additional cut, 9 e > 
8°(case-b). 
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Figure 8: The invariant mass distribution of q and q' calculated at yjs = Mz- 
The integrated luminosity and experimental cuts described in ref. [11] but 
applied at the generator level is assumed. 
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Figure 9: The transverse momentum of qq' system calculated at \/s — Mz- 
The integrated luminosity and experimental cuts described in ref. [11] but 
applied at the generator level are assumed. 



20 



25 E- 
30 



_L 



1 



Jj 







50 100 
M ud (GeV) 



150 



